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The mass and decay constant of the ground state of spin singlet bb is calculated 
by using the QCD sum rules method. This meson is an axial vector P-wave named 
hb(lP) meson with J p = 1 + . It is found that the mass of this meson is (m = 
9940 ± 37) MeV, which is consistent with the recent experimental data by Belle 
collaboration. 
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I. INTRODUCTION 

Heavy quarkonium bound state and their spectroscopy such as bb and cc have been 
considered as a good tool for the investigation of the Quantum Chromodynamics(QCD) 
interaction. Bottomonium is specially important because of the variety of bound states 
and decay channels allowing us to study and fix many parameters and aspects of Standard 
Model(SM) and QCD. Moreover, measurements of bottomonium masses, total widths and 
transition rates serve as important benchmarks for the predictions of QCD-inspired potential 
models, non-relativistic QCD, lattice QCD and QCD sum rules [l|, Hj. The spin-singlet P- 
wave bound states of bb such as hb{lP) can be used as tests of the P-wave spin-spin (or 
hyperfine) interaction. Therefore, theoretical calculations on the physical parameters of this 
meson and their comparison with experimental data could give essential information about 
their nature and hyperfine interaction. 

The QCD sum rules approach which enjoys two peculiar properties namely its 

foundation based on QCD Lagrangian and free of model dependent parameters, as one 
of the most well established non-perturbative methods has widely been applied to hadron 
physics. In the present work, we extend the application of this method to calculate the 
masses and decay constants of the axial vector P-wave hb(lP) meson with the quantum 
numbers J p = 1 + . Since the heavy quark condensates are suppressed by the inverse powers 
of the heavy quark mass [§, Toj | . so as the first nonperturbative contributions, we take into 



account and calculate the two-gluon condensate diagrams. Note that, Belle collaboration has 
recently reported the observation of the h b {lP) spin-singlet bottomonium state produced in 
the reaction e + e~ — > hb(lP)n + 7r with significances of 5.5a . They found that m[hb(lP)] = 
(9898.25 ± 1.061};$) MeV Q. 

The layout of the paper is as follows: in the next section, the sum rules for the mass and 
decay constant of the ground state axial vector meson is calculated. Section III encompasses 
our numerical predictions for the mass and leptonic decay constant of the hb(lP) axial vector 
meson. 



II. QCD SUM RULES FOR THE MASS AND DECAY CONSTANT OF THE 

h b (lP) AXIAL VECTOR MESON 

The two point correlation function respunsible for the mass and decay constant of the 
axial vector meson can be written as: 

K = if d 4 xe^(0 | T (#(*)#(())) I 0), (1) 

where T is the time ordering product and = b{x)^^^b{x) is responsible for creating the 
axial vector quarkonia from the vacuum with the same quantum numbers as the interpolating 
currents. 

According to the QCD sum rules, the correlation function has to be calculated in two 
different ways. Regarding the physical or phenomenological side, it can be evaluated in 
terms of hadronic parameters such as meson mass and decay constant. On the other hand, 
in the QCD or theoretical side, it is calculated in terms of QCD degrees of freedom such 
as quark masses and gluon condensates by means of the operator product expansion (OPE) 
in deep Euclidean region. Matching these two representations of the correlation function 
via the dispersion relation, the relation among the physical or hadronic parameters i.e., the 
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mass and the decay constant and the fundamental QCD parameters can be achieved. The 
Borel transformation with respect to the momentum squared is applied for both sides of the 
correlation function. This end is necessary to suppress the contributions of higher states 
and continuum. 

Firstly, to calculate the physical part a complete set of intermediate states with the same 
quantum numbers as the interpolating current is added. Performing the integral over x and 
isolating the ground state it is obtained that: 



n 



(0\JA(P)\A)(A\J?(P)\0) 



+ •••, (2) 

m A — p* 

where • • • represents contribution of the higher states and continuum and m A is the mass of 
the axial vector meson. 

The matrix element, parametrized in terms of the leptonic decay constant, is as follows: 



(0 | J(0) | A) = f A m A e ul 



(3) 



where f A is the leptonic decay constant of the axial vector meson. The summation over the 
polarization vectors are: 



-9^v + 



P»Pv 



771 



2 ' 



(4) 



Using Eq. (j3J), the final expressions of the physical side of the correlation function is as: 



n 



ft i' 



p z 



PjJK 
m 2 A 



(5) 



to proceed with the calculation of the mass and the decay constant the g uu structure is used. 







FIG. 1. (a): Bare loop diagram (b): Diagrams corresponding to gluon condensates. 



In QCD side, the perturbative (short distance) and non-perturbative (long distance) 
parts of the correlation function must be separated. In this regard, the OPE method is used 
and the correlation function is calculated in deep Euclidean region, p 2 <C —Aq CD . So, the 
correlation function is written as: 



u QCD = u pert + n 



'-nonpert 



(6) 



4 



The short distance contribution (bare loop diagram in figure flTJ part (a)) is calculated using 
the perturbation theory, whereas the long distance contributions (diagrams shown in figure 
(0Q) part (b) ) are parametrized in terms of gluon condensates. To proceed, we write the 
perturbative part in terms of a dispersion integral as follows: 

n J — ^2 ^-nonpert j (7) 

where, p(s) is the spectral density. This spectral density can be calculated having the 
Feynman amplitude of the bare loop diagram and using the Cutkosky rules. Note that, 
in the Cutkosky rules the quark propagators are replaced by Dirac delta function, i.e., 
_±_^ { _2m)5(p 2 -m 2 ). 

Finally, we find the spectral density as: 

The non-perturbative part is calculated via the gluon condensate diagrams represented in 
part (b) of figure ([T]). The Fock-Schwinger gauge, x^A a ^{x) = 0, is chosen. Then, the vacuum 
gluon field is as: 

Al(k') = - l -(27r)*G^(0)^\k'), (9) 
where k' is the gluon momentum and the quark-gluon-quark vertex is used as: 

r ^ = ^(y) , (10) 
As a result, the non-perturbative part in momentum space is obtained as: 



rl @ 
Konpert = J Q («^ 2 ) _ ^ - dx. (11) 

The explicit expressions for is as follows: 

6= ~x 2 ^2mtx 3 (m 2 b (4-3x) +p 2 (x 2 + x - 2)) 

+mtx(m 2 b (x(n - 2x(18(x - 3)x + 47)) + 8) + p 2 x(x(21x - 25) - 7)(x - l) 2 ) 
+m 2 b (-m 2 b p 2 (x - l)x(x(x(2x(27x - 82) + 149) - 32) - 11) 
+m 4 b (3x(x(x(3x(Ax - 17) + 76) - 46) + 8) + 5) + p\x - 1) V(24x 2 - 22x - 5)) 
+m 2 (x — l)x 2 (m 2 p 2 (7 — x{2x + 3)) + m b (3x — 5) 

+p\x - l)(2(x - l)x + 3)) + {x- l){-m 2 b p\x - l)x(2x(x(2a:(6x - 19) 
+37) - 10) - 3) + mlp 2 (x{x(x(x{27x - 112) + 162) - 97) + 14) + 5) 
+m£(x(57 - 2x(3x(2x - 9) + 43)) - 15) + p 6 (x - 1) V(6(x - l)x - 1)) 

+3mlx A + 3ml(x - 1) V(4x +1)1. (12) 
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The next step is to apply Borel transformation over p 2 for suppressing the contribution 
of the higher states and continuum. Then, sum rules are obtained as: 



m\r A e " 



4ml 

b 



ds p(s) e m 7 + BTl 



nonpert j 



(13) 



where M 2 and So are the Borel mass parameter and the continuum threshold, respectively. 

The meson mass (itia) is extracted with simple algebra (applying derivative with respect 
to — jf2 to the both sides of the above sum rule and dividing by itself). It is obtained as 
follows: 



m , 



JZlds P A (S) e~Th + BYltnpert 



£2 ds p A (s) e + BYi A onpert 



(14) 



where 



BYi 



nonpert 



A 



1 rn b 

g M 2 x(x — 1) 

n 7r96M 6 (x 



1)% 3 



(a s G 2 )dx, 



and we get 



A= -m\{x - l)x 2 (m 2 b (2x 2 - 5x + 2) 



(15) 



(16) 



+2M 2 x(x 2 - Ax + 3)) - m 4 b {x - l)x 3 (m 2 b (3x 2 - Qx + 2) 



llx 2 + 17a; - 6) 



+M 2 x(3x 2 -Ax + 1)) + M 2 x(2x ...... .... 

+ml(x 3 - Ax 2 + Ax - 1) + 2M 4 x 2 (2x 3 - Ax 2 + 5x - 3)) + m 2 b x(rn 2 b M 2 x{<ox 



-26a; 4 + 43a; 3 - 31x 2 + 9x - 1) + m 4 (3x 5 - 15x 4 + 27x 



+4MV(3x + l)(x - l) 4 ) + (x 



20x 2 + 7x - 1) 
l)(-2m 2 b M 4 x 3 (6x 4 - 28x 3 + A5x 2 - 31x + 8) 



+m 4 M 2 x(-3x 5 + 16x 4 - 33a; 3 + 28x 2 - 11a; + 2) - mt(x 5 - 6a; 4 + 13a; 3 - llx 2 
+5x - 1) + 2Af 6 (a; - l) 3 a; 3 (6a; 2 - 6a; - 1)) + m 6 b (x - l)x 5 + m 6 b (x - l) 2 x 3 . 



III. NUMERICAL RESULTS 

The sum rules for the decay constant and the mass of the hj,(lP) depend on the bottom 
quark mass, QCD parameters and two auxiliary parameters namely Borel mass parame- 
ter M 2 and continuum threshold s . The values m b = 4.8 GeV and (0 | \a s G 2 \ 0) = 
0.012 GeV 4 [ll[ are chosen. Even though the physical observables suppose to be inde- 
pendent of auxiliary parameters, the weak dependency of physical observables on auxiliary 
parameters can be tolerated. We have analysed the dependence of the mass and the decay 
constant on the Borel mass parameter M 2 and the continuum threshold so as shown by figs. 
(2) and (3). With these figures, the interval of auxiliary parameters can be fixed in principle 
when the sensitivity to the variation of these parameters are weak. Furthermore, the interval 
of the Borel mass parameter is determined when the higher state and continuum contribu- 
tions and the contributions of the highest order operators are suppressed; i.e., the sum rules 
for the mass and the decay constant converge. As a result of the procedure mentioned above, 
the interval of the Borel parameter is taken to be 15 GeV 2 < M 2 < 30 GeV 2 . Note that the 



similar region has been found for analysing the mass and the decay constant for A^o [12 



The continuum threshold So is not completely arbitrary but it is related to the energy of the 
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first exited states (hb(2P)) with the same quantum numbers as the interpolating currents. In 
other world, the interval of Sq must be in between the mass square of (/i&(1P)) and (hb(2P)), 
i.e., m 2 [h b (\P)} < s < m 2 [h b (2P)], where m[h b (lP)} = (9898.25 ± 1.06tJ$) Me V and 
m[h b {2P)} = (10259.76 ± 1.064t}$j) MeV Then, the interval 101 GeV 2 < s < 103 GeV 2 
for the continuum threshold So is chosen. Here, it is worth to be mentioned that in the 
above regions for Borel mass parameter M 2 and continuum threshold so the dependence of 
the results on auxiliary parameters is the most weak one. 

Regarding the intervals of the Borel mass parameter M 2 , continuum threshold so and 
input parameters mentioned above, the following values for the mass and the decay constant 
of /ib (IP) are obtained: 

m[h b (lP)}= (9940 ± 37) MeV (17) 
f[h b (lP)}= (94 ± 10) MeV (18) 

Note that, we use the experimental value of the mass m[h b (lP)] = 9898.25 MeV rather than 
the theoretical value of that in order to determine the decay constant. It is also worth to 
mention that the result for the value of the mass is in good agreement with the recent data 
on the mass of h b (lP) by Belle collaboration [l[ where they found: 

m[h b {lP)) = (9898.25 ± 1.06^) MeV (19) 

The mass and the decay constant of axial vector h b (lP) are depicted in figures 0U|5]) at 
three fixed values of the continuum threshold in terms of Borel mass parameter M 2 . These 
figures indicate a good stability of the mass and decay constant with respect to the Borel 
mass parameter M 2 . From these figures, it is also clear that the results depend weakly on 
the continuum threshold in its working region. 

In summary, this study presents first theoretical calculation on the mass and decay con- 
stant of the newly found axial vector meson h b (lP) in the framework of the two point QCD 
sum rules. While the result for the mass is consistent with the measured value by Belle 
collaboration, the decay constant has not yet been measured. 



I. Adachi et at [Belle Collaboration]. l arXiv:1103.3419l [hep-ex]. 
N. Brambilla et al, Eur.Phys.J.C 71, 1534 (2011). 

M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, Nucl. Phys. B 147 (1979) 385. 

P. Colangelo, A. Khodjamirian, in "At the Frontier of Particle Physics/Handbook of QCD", 

edited by M. Shifman (World Scientific, Singapore, 2001), Vol. 3, p. 1495. 



V. M. Braun, arXiv: hep-ph/9801222 



I. I. Balitsky, V. M. Braun, A. V. Kolesnichenko, Nucl. Phys. B 312 (1989) 509. 
T. M. Aliev, K. Azizi, V. Bashiry, J. Phys. G 37, 025001 (2010). 

V. Bashiry, K. Azizi and S. Sultansoy, Phys. Rev. D 84, 036006 (2011) [ arXiv: 1104.28791 
[hep-ph]]. 

S. Narison, Phys. Lett. B387, 162 (1996). 

M. Jamin, A. Pich, Nucl. Phys. B507, 334-352 (1997). [hep-ph/9702276] . 
B. L. Ioffe, K. N. Zyablyuk, Eur. Phys. J. C27, 229-241 (2003). [arXiv:hep-ph/0207183|| [hep- 

Ph]]- 

[12] E. V. Veliev, H. Sundu, K. Azizi, M. Bayar, Phys. Rev. D82, 056012 (2010). [ arXiv:1003.0lT9l 
[hep-ph]]. 



7 




FIG. 2. Dependence of the mass of the axial vector 67^756 on the Borel parameter, M 2 and the 
continuum threshold sq. . 




FIG. 3. Dependence of decay constant of the axial vector 67^756 on the Borel parameter, M 2 and 
the continuum threshold sq. . 
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FIG. 4. Dependence of the mass of the axial vector 67^756 on the Borel parameter, M 2 at three 
fixed values of the continuum threshold. The upper, middle and lower lines belong to the values 
103 GeV 2 , s = 102 GeV 2 and s = 101 GeV 2 , respectively 
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FIG. 5. Dependence of the decay constant of the axial vector bb on the Borel parameter, M 2 at 
three fixed values of the continuum threshold. The upper, middle and lower lines belong to the 
values 103 GeV 2 , s = 102 GeV 2 and s = 101 GeV 2 , respectively. 



